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Conventional methods to determine non-perturbative parameters in QCD, such as the different 
variants of QCD sum rules or the minimal hadronic approximation, combine a certain degree of 
matching to QCD with inputs from hadronic parameters. The latter introduce systematic errors 
difficult to quantify. In this paper I will apply a method based on rational approximant theory 
where matching is maximized and no hadronic inputs are used, thereby leading to simple analytical 
relations between high and low energy parameters. I will be mostly interested in the phenomeno- 
logical applications to the Ulr and Hvt correlators, with especial emphasis on quantities like the 
d = 6 and d — 8 vacuum condensates in Ulr or the quark condensate magnetic susceptibility \o- 



INTRODUCTION 



Confining gauge theories like QCD are extremely 
challenging. Information based on first principles is 
only available at high space-like momenta, where per- 
turbation theory is valid, while at low energies the 
spontaneous breaking of chiral symmetry allows a de- 
scription in terms of an effective theory of pions known 
as chiral perturbation theory (ChPT). 

In order to go beyond these two energy regimes, 
in the so-called resonance region, new genuinely non- 
perturbative techniques have to be employed. In this 
respect, lattice QCD is probably the best-positioned 
candidate to eventually unravel non-perturbative 
physics, but even if this goal is accomplished, it should 
be complemented with analytical non-perturbative 
methods. 

Over the years many different methods have ap- 
peared to estimate non-perturbative parameters in 
QCD. Among them, the most popular are perhaps 
vector meson dominance (VMD) [lj, QCD sum rules 
(QCDSR) with all its variants - finite-energy sum 
rules (FESR), Borel sum rules (BSR), etc.- or the min- 
imal hadronic approximation (MHA) 0, Q • 

Despite the differences between those approaches, 
there is a common ground: they all adopt a mero- 
morphic hadronic ansatz, i.e. some modeling of the 
hadronic spectrum consisting of single poles. This can 
be connected to a well-defined limit of QCD, namely 
the limit of large number of colors Q. In this limit, 
QCD looks like a topological field theory which admits 
a power expansion in the parameter 1/N C . At lead- 
ing order, any correlator is a meromorphic function. 
Plenty of phenomenological features of QCD can be 
qualitatively understood by inspecting the leading or- 
der in the l/N c expansion [6(. However, a quantitative 
solution to large- N c QCD has proved elusive so far. In 
other words, the values for the (infinite) poles and their 



associated residues remain unknown. 1 

Therefore, from this perspective, VMD, QCDSR and 
the MHA may be seen as different phenomenological 
approximations to QCDoo. VMD assumes that, when- 
ever vectors contribute to a certain process, the first 
vector resonance in the spectrum gives the bulk of the 
non-perturbative effects. A natural extension of VMD 
is lowest meson dominance (LMD), which states that 
each channel is dominated by the lowest lying reso- 
nance. The problem is that very little is known about 
poles and residues of resonances. QCDSR, in contrast, 
is a method that allows to determine resonance param- 
eters by matching the hadronic ansatz to the operator 
product expansion (OPE). The MHA originated from 
a completely different context, namely the computa- 
tion of electroweak observables, like the electromag- 
netic pion mass difference Q, which can be expressed 
in general as integrals over the Euclidean regime of 
QCD correlators. The MHA also relies on a matching, 
but one typically matches to high and low energies. 
Therefore, VMD has no matching, QCSR is a 1-point 
matching procedure and the MHA a 2-point one. 

Only recently it has been pointed out [9( that the 
previous methods are special cases of what is known 
in the mathematical community as rational approxi- 
mants. Essentially every method based on meromor- 
phic ansatze and local matching is a Pade approximant: 
a Pade-type if the poles are fixed to physical masses, a 
partial-Pade if a subset of poles and residues is taken 
from experiment, or a plain Pade approximant if no in- 
put other than matching is used. Therefore, the MHA 
is a 2-point Pade-type approximant, FESR are I-point 
Pade-type approximants while BSR are 1-point Borel- 
Pade approximants. Only VMD falls out of this cate- 
gory because there is no matching involved. 



1 See however Ref. 0] for a strategy to determine a set of decay 
couplings in the large-7V c limit. 
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There is one important remark to be made: the 
theory of Pade approximants states clearly that the 
parameters that come out of matching (poles and 
residues) are not the physical ones, and only in cer- 
tain cases can be identified as such. In fact, as we will 
see later on, it is not uncommmon that some of the 
mass-parameters become complex-valued. As a result, 
a naive identification of Pade parameters with physical 
ones should be avoided. 

An additional and closely related characteristic of 
Pade approximants is that Euclidean quantities are 
quite insensitive to the details of the hadronic ansatz. 
This was already observed in Ref. [l(| HH in the con- 
text of MHA, where it was shown that meromorphic 
correlators, which bear little resemblance to QCD in 
the physical axis, can nonetheless be made extremely 
accurate in the Euclidean regime, provided the corre- 
lator complies with known properties of QCD at high 
and low energies. This was a direct consequence of the 
MHA being a Pade approximant, though at the time it 
was not fully realized. Therefore, if one is interested in 
Euclidean parameters there is no need to use hadronic 
input. 

In this paper I will explore the consequences of work- 
ing with a minimal (meromorphic) ansatz for correla- 
tors constrained in a maximal way. This maximally- 
constrained approach differs from conventional MHA 
in that not only the decay couplings but also the masses 
will be left free. Thus, the method is none other than 
a plain Pade approximant: all parameters of the in- 
terpolator will be determined using only information 
from Euclidean space, with no reference to hadronic 
parameters whatsoever. In this work I will concen- 
trate mostly on the two-point correlators Ulr and 
ILvt- The Pade approximant will lead to definite rela- 
tions between low and high energy parameters, allow- 
ing in particular to relate the ChPT parameters Liq 
and Cg7 to the dimension-six and dimension-eight vac- 
uum condensates, or the magnetic susceptibility %o to 
the mixed and quark condensates. 

The structure of the paper will be as follows: in Sec- 
tion [H] I will discuss the application of the method to 
the Ulr correlator and its implications for the d = 6 
and d = 8 vacuum condensates. In Section [Ml I make 
contact with the theory of Pade approximants and 
show how the method is embedded in it. In Section HVl 
I extend the analysis to the correlator to give a 
prediction of the low energy parameter xo, the mag- 
netic susceptibility of the quark condensate, in terms 
of the mixed condensate. In Section [V] I briefly com- 
ment on the implications of this strategy for the scalar 
and tensor two-point functions. Conclusions and fu- 
ture prospects are summarized in Section IVTl Finally, 
I include an Appendix with the d = 6 and d = 8 four- 
quark operators for the vector, scalar and tensor sec- 
tors together with their expressions in the factorization 
approximation. 



II. OPE CONDENSATES IN U LR 

Consider the following two-point correlators 

n£T(?) = i J d A xe^(0\T{V^x)V^0)}\0) , 

K A (l) = * J d 4 ze l ^(0|T{A^)4(0)}|0) , 

(1) 

where V^{x) = u(x)^/^d(x) and A^(x) = u{x)^^d{x). 
We will define its difference as Ulr : = Hvv — TIaa- In 
the chiral limit, Lorentz and gauge invariance lead to 

n^(g) = W - q 2 gn^LR(q 2 ) ■ (2) 

This correlator has a lot of interesting properties. First 
of all, it is an order parameter of the spontaneous 
breakdown of chiral symmetry (SxSB). In particular, 
this means that it vanishes to all orders in perturbation 
theory. Therefore, all its high and low-energy coeffi- 
cients are order parameters of SxSB too. Second, its 
first two OPE terms cancel on general grounds: gauge 
invariance forbids a, d — 2 condensate, while the d = 4 
operator is purely gluonic (in the chiral limit) and 
therefore chirally symmetric. This led to the celebrated 
superconvergence relations between hadronic parame- 
ters for the vector and axial channels [l2T |. nowadays 
known as Weinberg sum rules. The high energy fall-off 
can therefore be cast as 

oo f 

H y ' n=3 H 

On the other hand, ChPT allows to parameterize the 
low energy regime as 

lim H LR (q 2 ) = 4 - 8£io + 16CW + 0(q 4 ) , (4) 

q 2 ^Q q A 

where Lio and Cs7 are low energy parameters, encod- 
ing the dynamics of the hadronic spectrum. 2 

When the number of colors is large, any correlator 
is saturated by the exchange of an infinite number of 
stable one-particle states. For Ulr, the absorptive part 
reads 

. oc 

-lmU L R(t) = -f^(t)+J2fvJ(t-m 2 Vn )- 

7T z — J 

n 

oo 

J2fln S (t-<n) , (5) 



2 In general, both the vaccum condensates and ChPT parame- 
ters denned in Eqs. J3jl and (J4j have a logarithmic dependence 
when quantum corrections are taken into account. In this pa- 
per I will always consider their leading-order scale-independent 
contribution. 
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where we have defined 



\V M \p n (p,X)) = fvnmvn^ 

I Aft\ p n (jp, A)) = fAnniAr 



(A) 



(6) 



Iliij satisfies an unsubtracted dispersion relation and 
therefore one can write 
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If the previous expression is regarded as an interpolator 
of the true QCD correlator, then the unknown masses 
and decay couplings can be determined, for instance, 
by matching to the known OPE expansion. In general, 
for an arbitrary large but finite spectrum, one finds the 
following set of matching equations 
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One of the main advantages of meromorphic approx- 
imants is that low and high energy parameters can be 
related analytically. Following our discussion in the In- 
troduction, I will restrict my attention to the minimal 
number of states per channel, 
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(10) 



In order to fully determine masses and decay constants, 
we need to solve the following system of equations: 



rZ z 

Ja^a 
pA m A 



/a — fv — 

4-2 2 

fv m v = 



j-2 4 

fWv 



, 

^8 , 



(11) 



which is a particular case of Eq. ([8]), where the first 
two equations are the celebrated Weinberg sum rules. 
With the masses free, the previous system of equations 
becomes non-linear, and its solution becomes consid- 
erably involved. However, it is instructive to inspect 
the solution at a qualitative level. The first thing to 
notice is that there are actually two solutions, as a con- 



sequence of the symmetry {f v 



of the ansatz. The second and most relevant fact is 
that the parameters can become complex if 
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M v 



2 j-2 
l Vn 



>2j 



(8) 



whose solution determines the hadronic parameters in 
terms of vacuum expectation values of QCD operators. 
The previous system of equations was studied in de- 
tail in Ref. [l3l |. where very interesting remarks were 
made on the relation between OPE condensates and 
the pattern of vector and axial states in the spectrum 
of Ulr- In particular, the structure of Eqs. ([5]) leads 
to a Vandermonde-type linear system and can there- 
fore be solved analytically. The results show that all 
OPE condensates, as defined in Eq. ([3]), should be posi- 
tive. This complies with a theorem by Witten [lj] that 
states the following positivity condition: 



q 2 K LR {q 2 ) > 0, -00 < q 2 < 



(9) 



As pointed out in Ref. [131 ] , Witten's inequality in con- 
junction with the Weinberg sum rules requires ^6 to be 
positive. 

Note nonetheless that the conclusions reached in 
Ref. 13] hold if the masses in Eq. ([5]) are taken to be 
physical, i.e., if one is constructing a Pade-type. How- 
ever, as I discussed above, plain Pade approximants 
do not impose such constraints and instead leave all 
masses and decay constants as free parameters. In the 
following, I will explore the consequences of analysing 
Eqs. © as a plain Pade. 



e 8 < 4/- 2 £ 6 3 



(12) 



Imaginary solutions obviously invalidate any interpre- 
tation of the parameters as being physical. However, 
notice that this does not affect Euclidean quantities. 
For instance, if we define the parameters of the chiral 
expansion as 



lim U LR (q 

qr 2 — >0 



= f+Ec^> as) 



the solution of Eqs. (fTTj) turns out to be amusingly 
simple and expressible as a general recursive relation: 
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(14) 



where Uj (x) are Chebyshev polynomials of the second 
kind. The first two chiral coefficients defined in Eq. (H} 
therefore adopt the compact expressions: 
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(Zip 



e 6 ) 



(16) 



Let us concentrate for now on the first relation, which 
is one of the main results of this paper. Notice in the 
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first place that, since Liq < 0, it predicts the sign of £s 
to be positive. This can now be compared to existing 
determinations of the condensates. 

In the first two columns of Table I I list the values 
for £ 6 and £s reported in Refs. [IB]- [HI- Note that 
Eq. ([T5|) is only compatible with the first half of the 
table, i.e., those determinations where £s > 0. The 
third column lists the values that £s would take, in the 
different analyses, if L\a and £g were taken as inputs 
according to Eq. (Ti"5|) . We adopt the most recent value 
for L\q [25| 

L w (m p ) = -(5.22 ±0.06) • 1(T 3 ; (17) 

extracted from a 0(p A ) analyses of tau decay data, 
together with f n = (130.4 ± 0.2) MeV. 

Notice the remarkable agreement between columns 
2 and 3 when Eq. (fT5|) is used. This is also illustrated 
in Fig. [TJ Therefore, Eq. (fl5|) not only predicts the 
sign for £s to be positive but seems to indicate that 
the first half of Table I can be described to a very good 
approximation by the quadratic relation 



& = ce 6 , c 



-8ii 



f 4 

71 



(18) 



where the constant C is nowadays known to the 1% 
level. 

If we now put together Eq. (fTS"]) and Eq. (fl6|) one 
can give a prediction for the condensates. Inverting 
the system one finds 



f 

■I 71 



16(4L 2 - C s7 f 2 ) ' 



£s = - 



L wf% 



32(4L 2 - C s7 f 2 ) 2 ■ 



(19) 
(20) 



L\q and Cs7 could in principle be taken from the recent 
0(p 6 ) analyses on tau decay data [25j |: 

L w (m p ) = -(4.06 + 0.39) ■ 10" 3 , 

C S 7(m p ) = (4.89 ± 0.19) • 10~ 3 GcW 2 . (21) 

However, there are some reasons not to proceed this 
way. First and foremost, the previous values for the 
chiral coefficients would lead to a negative dimension- 
six condensate in Eq. (fTl?)) , contradicting Witten's in- 
equality. This can be easily seen from the denominator 
of Eq. (frnj) • The positivity condition is 



087 < -72~ 

J 7V 



(22) 



Of course this does not mean that the values of Ref. [25| 
violate Witten's theorem. What it means is that there 
is an apparent incompatibility between the method 
developed in this paper and the results reported in 
Ref. [25|] . The origin of the discrepancy is unclear to 



me at this point. However, a bigger value for Liq in 
Eq. (f2~Tj) . closer to the 0(p 4 ) result would: (a) comply 
with the positivity condition and (b) solve an appar- 
ent puzzle: the 0(p 4 ) and 0(p 6 ) values in Ref. [25| are 
not compatible within errors. It is also significative 
that if Lio = —5.22 ■ 10~ 3 is used, the values of the 
condensates for the first half of Table 1 come out from 
Eqs. p^]) and (|2TJ)) with natural values for Cs7, in the 
range 0.0037 GeV~ 2 < C 87 < 0.0047 GeV~ 2 . 

An alternative way to determine the vacuum con- 
densates is through a fit with the values of the first 
half of Table I subject to Eqs. (HU) and OH). To be 
more precise, I will find the values for and £s that 
best fit Eq. (|16p while keeping Eq. (|15p as a constraint 
with L w given by Eq. ([17) . 3 

The results are 

£ 6 = (+7.6 ±0.4) • 10~ 3 GeV 6 ; (23) 
£ 8 = (+8.3 ± 1.0) • 10~ 3 GeV 8 ; (24) 

leading to a very reasonable value for Cs7- 

C s7 = (+4.0 + 0.2) • 10- 3 GcV~ 2 . (25) 

As a consistency check, one could now compute the 
prediction for the electromagnetic pion mass difference, 
which is related to Hlr as [8j 

Am - = -7^2 / dQ 2 Q 2 U LR (Q 2 ) . (26) 

With ansatze like Eq. (|10| satisfying the Weinberg sum 
rules, the previous equation can be expressed as 

. 3a m, 2 A my ( m 2 A \ 

Am, = / V 2 log -T ■ (27) 

87rm,r m A — my \my / 

Plugging in the values found in Eqs. and (|2"4"|) one 

gets 



Am, = (4.60 ± 0.28)MeV 



(28) 



in excellent agreement with the experimental value 
Am, = (4.5936 + 0.0005) MeV. 

I would like to point out that the values found above 
only include statistical uncertainties. Systematics both 
due to quantum corrections and intrinsic to the nature 
of the rational approximation as an iterative method 
are expected to be dominant. Their impact on the 
different parameters is difficult to estimate, but gener- 
ically one should expect a 10-30% correction. 



3 In other words, I assume that Lio is known with certainty. If 
both Eqs. (1 1 5 I t and I I16H were weighted evenly in the fit, the 
resulting value for Lio would have unrealistically large errors. 
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Cirigliano al. on OPAL [211 


+5.43 ± 0.76 


-1.35 + 3.47 


+4.3 + 1.2 


Bijnens et al. on ALEPH [22] 


+3.4t« 


-UA±l° * 


+1.7 + 2.4 


Bijnens et al. on OPAL [22] 


+4.0 ± 2.0 


-10.4l 8 6 ;° 
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TABLE I: Values for the dimension-six and dimension-eight OPE condensates (in 1CP 3 GeV 6 and 10 -3 GeV 8 , 

respectively) reported using different phenomcnological techniques. In the last column we list the would-be 
value for the dimension-eight condensate if Eq. (|15[) were used, taking as input the values for L 10 and /Z listed 

in the main text and the values for £ 6 from the first column. 



III. THE METHOD AS A PADE 
APPROXIMANT 



At this point, let me come back to an issue I have 
mentioned but not discussed. So far all the quantities 
evaluated, namely chiral coefficients, vacuum conden- 
sates and even the pion electromagnetic mass differ- 
ence, have been real, as one expects, yet the hadronic 
parameters from the ansatz, that one would naively 
identify with physical states, might become complex if 
the reality condition £g > 4/~ 2 £| is not fulfilled. In- 
deed, most of the values reported in the first half of 
Table I do not satisfy the reality condition, and yet 
they still comply with Eq. (|15p . Euclidean quantities 
seem to be mysteriously protected. 

This is where the theory of Pade approximants be- 
comes useful. As I already mentioned in the intro- 
duction, almost the entire analytical techniques that 
deal with non-perturbative strong interactions can be 
formulated as different kinds of Pade approximants to 
meromorphic functions. A Pade approximant is 
the ratio between two polynomials, where the indices 
n, m denote the degree of numerator and denominator, 
respectively. 

By definition, the approximant should be regular 
at the origin. For Wlr regularity at the origin can 
be achieved if one works with the function q 2 IlLR in- 
stead. This is a P% approximant which should be sup- 
plemented with four constraints, which for convenience 



will be rewritten as 

fv ~ fA = fir ) 
fv m V - fA m A = , 

fv Ia _ or 
2 2~ — _ O-Lio , 

m v m A 

= 16C 87 . (29) 

rriy m A 

The effect of the two Weinberg sum rules is to reduce 
the approximant to a P§ approximant of the form 

„2tt _ m V m Afl /o n ^ 

q H* 2 r^) ■ (30) 

In a regular Pade-type approximant the masses in the 
previous formula can be taken from experiment. In 
contrast, in a plain Pade approximant two more con- 
straints are needed such that the masses are functions 
of Lio and C%-j. 

Pade approximants are iterative methods, and the 
exercise done in the previous section is just the first 
iterative step. In order to improve the approximation 
there is a natural and well-defined prescription: new 
resonance contributions can be added to the ansatz, 
and for each one two more constraints from low ener- 
gies should be added to Eqs. ([2^|) . 4 The iterative proce- 



4 Given that there is precise data on the absorptive part of II 
this is actually feasible and will be the subject of a separate 
article. 
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dure thus defined is empowered by the following result: 
when the function is meromorphic, there is a theorem 
by Pommerenke [26[ that ensures convergence of as 
n, m — > oo on a compact subset of the complex plane 
except on a set of null measure. 5 This includes the 
physical axis, but also spurious poles that the Pade 
approximant can generate. These spurious poles can 
be imaginary, in which case they always come in com- 
plex conjugate pairs. 

For instance, for the values of Eqs. ([23]) and ([24]) . one 



finds m% = (m 2 A )* = (0.546 + 0.386i)GeV 2 , in agree- 
ment with Pade theory. There is nothing pathological 
about the appearance of complex poles: it is the price 
to pay for approximating an infinite number of poles 
by a finite number of them. One should keep in mind 
that the parameters in the ansatz are not physical and 
should be seen instead as the effective masses and de- 
cay couplings that best approximate the function as a 
whole. As an illustration of this last point, consider 
the pion electromagnetic mass difference in Eq. (|2"T|) . 
Using the MHA, rriA, my should be identified with the 
physical masses for the first vector and axial hadronic 
states. This would yield Am, ~ 6.0 MeV, while the 
maximally-constrained approach yields Am, ~ 4.60 
MeV. Despite the obvious differences in the interpre- 
tation and values for the mass parameters, in this case 
both methods yield a consistent prediction for Am,. I 
cannot think of a better example to illustrate the na- 
ture of Pade approximants and the fact that radically 
different ansatze in Minkowski space give reliable pre- 
dictions in Euclidean space. 

This last example brings up an interesting question, 
namely how our the results of the last section would 
change if different Pade approximants, like Pade- type 
or partial-Pade approximants, were used. For in- 
stance, one could add a new term in the ansatz of 
Eq. (fT0|) . while fixing the first two poles to the p(770) 
and ai(1260) physical masses, and solve the system 
of constraint equations again (partial-Pade). Or, al- 
ternatively, one could apply an MHA-like ansatz by 
adding four resonances, fixing their poles to the first 
four experimental masses in the Hl^ spectrum, and 
use the constraints to solve for their decay couplings. 
It turns out that in the partial-Pade case the multiple 
solutions for the matching equations, unlike the plain 
Pade approximant, do not lead to unique predictions 
for ^6 or £§. In contrast, the use of a MHA-like ansatz 
leads to a very stable Lio, even if the values for the 
masses are drastically changed. However, the values of 
the vacuum condensates are extremely sensitive to such 



changes. This sensitivity of the MHA to high energy 
parameters was already observed in Ref . Q . 



IV. THE MAGNETIC SUSCEPTIBILITY OF 
THE QUARK CONDENSATE 

We can apply the strategy described in the last sec- 
tions to obtain a prediction for the magnetic suscep- 
tibility of the quark condensate xo, a quantity intro- 
duced in Ref. [27| in the study of the scattering of nucle- 
ons in electromagnetic backgrounds. Since then, xo na s 
been subjected to multiple determinations, but there 
is still controversy as to its value, with many determi- 
nations hovering around xo ~ 3 GeV~ 2 , while others 
pointing at higher values, xo ~ 9 GeV -2 . Recently, it 
has been shown that the quantity plays a role in the 
hadronic lig ht-by-light scattering contribution to the 



(5-2), 



[29j , adding an extra motivation to resolve 



the present discrepancy. 

I start considering the following two-point function 

n^ p (g) = l U A x e**{0 | T{ V,{x) 13,(0) }| ) , 



(31) 

where T vp {x) = u(x)a vp d(x), and a^ v = i/2 [7", 7"]. 
It is easy to check that Hvt is an order parameter of 
S^SB, which is non- vanishing because at the hadronic 
level fy n , defined as 

(0|T^| P „(p,A)) = iti n (eMpu - eWp p ) , (32) 

is non-zero. Using Lorentz, parity and gauge invariancc 
one can factor out the tensor structure as follows: 

n£7(?) - * {<f<T - 1 V 9 W ) Tlvriq 2 ) ■ (33) 
At high energies, the OPE expansion reads [3C||: 6 

lim U VT (q = 2— — + • • • , 34 

q 2 ^(-oo) q A 3 q 1 

where we used the short-hand notation G = G^a^ . 
At low energies the xo parameter is defined as 



lim U VT (q 2 ) = -Xo{tpi>) + 



(35) 



In the large- Nc limit, the absorptive part of the corre- 
lator takes the form 

^ 00 

- Imn yT (£) = V" fvnfvnm Vn 5{t - m Vn ) , (36) 



Notice that strictly speaking Pommerenke's theorem does not 
apply when q 2 —* 00, i.e., when dealing with vacuum conden- 
sates. However, numerical exercises with toy models seem to 
find convergence also there. We refer the reader to Ref. [t| for 
details. 



6 Note that our sign convention for g s differs from the one em- 
ployed in [30j . 
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FIG. 1: Values for the dimension-6 (in GeV 6 ) and dimension-8 (in GeV 8 ) condensates from the first two 
columns of Table I (only points with positive £s have been included), and its comparison with Eq. (fTS")) . 



and therefore the correlator can be expressed as 

fvn m Vn 



n 



VT 



- m 



Vn 



x _ fvn 
An — — 

JVn 



(37) 



which follows again from the fact that ILvt satifies an 
unsubtracted dispersion relation. 

With the information given in Eq. (|34|) . our ansatz 
for the spectral function will consist of a single reso- 
nance state. Matching the high-g 2 limit of Eq. (f3"7| to 
the OPE, the resulting constraint equations are 



-2(#) , 

2g s 



(38) 



This system of equations is easily solvable, and it pre- 
dicts the following value for the magnetic susceptibil- 
ity: 



Xo = - 



(ipip)m v gsitpGip) 
It is customary to define the parameter 



mP. as 



(39) 



(40) 



whose value has been estimated from sum rules to be 
ml = 0.8 ± 0.2 GeV 2 [H- Subsequent lattice [U and 
sum rule analyses [l7T | are in very good agreement with 
this initial estimate, while instanton vacuum model cal- 
culations [H, [33[ predict slightly bigger results. 

Using the value for m\ reported in Ref. [3(3], we 
therefore find the prediction 



G 



Xo 



(7.5 ± 1.9) GeV" 



(41) 



which is big compared to the values obtained from 
QCD sum rule analyses (34j, instanton vaccum model 
calculations [35[ or exclusive B meson decays [3(1, but 
compatible with the values reported in Refs. [271 ] and 

The latter analysis is based on considerations on the 
transverse part of anomalous three-point functions and 
provides an analytical expression for xo- In our nota- 
tion, his result reads 



Nr. 



Xo 



27T 2 / 2 



8.9 GeV" 



(42) 



Combining the previous equation with the result in 
Eq. pij) . we can get an analytical prediction for the 
value of m 2 ,, namely 

m 2 Q = 4tt 2 / 2 ~ 0.67 GeV 2 , (43) 

which is entirely consistent with previous determina- 
tions. 

An additional non-trivial relation can be established 
between L\o and ml, thus checking the consistency 
of our approach between Ulr and Hvt, by express- 
ing Eq. (Ti"5")) under the factorization hypothesis. Using 
Eqs. (|A15[) and (|A16[) in the Appendix one finds 



Liq — — - 



)127ra, 



(44) 



The combination ira s {qq) 2 can be determined from the 
value of £ 6 from Eq. (f2"3"|) . The result is ■na s {qq) 2 = 
(10.7±0.6) x 10~ 4 GeV 6 , where the error is only statis- 
tical. This is in good agreement with the experimental 
value Tia s (qq) 2 ~ (9 ± 2) x 10~ 4 GeV 6 , meaning that 
our determination of ^6 is consistent with factorization. 

One could now solve for TOq in Eq. (|4"4"1) . With the 
experimental inputs for L\q and ira s (qq) 2 we find 

ml = (0.92 ±0.21) GeV 2 . (45) 
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This value can be compared with 



Scalar sector 



nin 



e 6 



f- = (1.13 ±0.06) GeV 2 



(46) 



showing a non-trivial global consistency between the 
experimental values for Liq and the quark condensate, 
our values for the vacuum condensates in Ulr and Hvt 
and the factorization hypothesis. 

Notice that according to Eq. (|4"Tj) . a smaller value for 
Xo requires a bigger value for the mixed condensate. 
For instance, in order to reproduce the sum rule value 
Xo = 3.15 GeV -2 [Hj], one would need ~ 1.9 GeV, 
certainly too big, even for instanton vacuum model es- 
timates. Our results therefore clearly favor a big value 
for xo- 

Finally, it is worth insisting that the parameter my 
has nothing to do with the my introduced in the pre- 
vious section when we discussed the vector channel. 
It is true that in the large-iV c limit the masses my n 
in Eqs. ([5|) and (l36|) are the same, but as soon as we 
truncate the spectrum the effective masses my and my 
become unrelated parameters. For instance, it is easy 
to verify that, while we saw that my was complex, my 
is real and given by 



my = —=mo ~ 516 MeV 
v3 



(47) 



Clearly, my should not be associated with the p(770) 
meson mass. Moreover, note that the following combi- 
nation of decay constants gives 



fyfv 



2V3( 



n) 



m 



(48) 



where we have used (go) = —(250 MeV) 3 and the val- 
ues reported in Ref. [38| for the decay constants. 



A COMMENT ON THE SCALAR AND 
TENSOR SECTORS 



The analyses we have performed in the previous sec- 
tions can be carried over to two-point functions involv- 
ing scalar, pseudoscalar and tensor currents. In the 
scalar sector, the lack of gauge invariance precludes a 
simple expression like Eq. (fTS"]) for the low energy cou- 
pling Lg and one is forced to resort to d = 10 vacuum 
condensates. The tensor sector, despite being out of ex- 
perimental reach, has attracted an increasing interest, 
especially in the study of nucleon structure. In partic- 
ular, the so-called tensor susceptibility has been 
estimated using sum rules and there is controversy as 



to its magnitude. According to some analyses [40l. I4l| . 
the discrepancies might be due to the hypothesis of 
vector meson dominance being invalid. 



Consider the following two-point correlators: 
liss(q) = * I d 4 xe^ x (0\T{S(x)SHo)}\0) , 



Upp(q) = i J d i xe^(Q\T{P(x)PHQ)}\0) , 

(49) 

where S^x) =: u(x)d(x) : and P^{x) —: u(x)ij5d(x) :. 
Analogously to the vector-axial case, we will be inter- 
ested in the difference ITs_p = IIss — IIpp, which is 
an order parameter of SxSB. Its operator product ex- 
pansion is of the form 



lim U S -p{q ) 

-►(-°o) 



n— 3 
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,2n— 2 



(50) 



while the low energies can be parameterized as follows: 

Bin 



lim Us-p(q 

q 2 ^0 



(51) 



Its spectral function in the large-iV c limit takes the 
form 



~Imn<j_p(t) = 

7T 



-B 2 f^(t)+J^fl n S(t-m% n ) 



fpj(t 



m 



Pn 



(52) 



and, since the correlator satisfies an unsubtracted dis- 
persion relation, one ends up with 



n 



S-P 



(<z 2 ) = 



r 



n 



Sn 



Pn 



m. 



In order to fulfill the minimal ansatz, we need four 
high-energy constraints. The matching equations take 
the form 



fp 



fs 



BlJl , 



j-2 2 

f P m P 

y2 4 

f P m P 
fpm 6 P 



r2 2 

fs m s 
fs m s 
f S m% 



do 



(54) 



Notice, in contrast to the vectorial case, that here we 
need to go as far as d — 10 operators. This is a con- 
sequence of the form of Eq. (fST)]) , which eventually can 
be traced back to gauge invariance: in the vectorial 
channel the Ward identity makes Ulr dimensionless, 
while n,g_p is a dimcnsion-2 object. One can solve 
Eqs. ([54"]) for the hadronic parameters in terms of OPE 
condensates to get the relation 



1 



32^2 



32B 2 



fp 



e 8 - tew 



(55) 
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With the previous equation, and taking as input the 
values for Lg, ^6 an d £s (assuming factorization), one 
can give a prediction for £10 • Using the results of 
the appendix and Lg — (9 ± 3) ■ 10~ 4 , one finds 
£10 = (0.002 - 0.02) GeV 10 . Recent results from the 
lattice UK seem to suggest much lower values for Lg. 
In particular, L%{m p ) = (3 — 4) • 10~ 4 . If this is con- 
firmed then our prediction for £10 would become nega- 
tive, £10 = -(0.002 - 0.003) GeV 10 . 

Again, it is worth stressing that the parameters on 
the left hand-side of Eqs. ([54")) are not to be taken as 
the masses and decay constants of physical particles. 7 
In fact, for typical values of the condensates and L% 
one always finds that at least one pole and one residue 
are complex. 



B. Tensor sector 



At low energies, the correlator takes the form 

n TT ( g 2 ) = 2^ + 2n 51 + = 2 Aj _ 32A , 

q z 2 q z 

_ (60) 

where we used the conventions of Ref . [46[ . The spec- 
tral function in the large- N c limit is given by 

1 x x 

- ImlT TT = ]T f 2 Bn S(t-m% n )-J2(fvn) 2 S(t-m 2 Vn ) , 

n n 

(61) 

where B n are (1 + ) mesons and fs n is defined as 

(0| V| B n (p, A)) = if Bn e^r, P eJ x) P p ■ (62) 

Knowing that Utt satisfies an unsubtracted disper- 
sion relation, the correlator with its minimal hadronic 
content reads 



Consider the two-point correlator 



n 



TT 



(q)=i d 4 «^(0|T{^(x)T^(0)}|0! 



(56) 

where, as before, T^^x) = u{x)a ilv d{x). Lorentz in- 
variance and the antisymmetry of tensor indices lead 
to the following decomposition 



IT 



TT 



(q) 



ii TT ( q 2 )F^ a(3 ( q ) + 



U+ T {q 2 )F^{q) , 
(57) 

where, for phenomenological purposes, we have con- 
veniently projected the form factors in combinations 
with well-defined parity, IIy T . F^' a ^ and F^ u ' a/3 are 
Lorentz tensors given by 



F 



fiu,af3 



<f<fa va 



(q) 

Ft 



Ft" (<il = -e'—'-e^-- g ni q P q T 

- ^ af \q) + q 2 {sTiff* - g^g va ) . (58) 



In fact, they are chiral projectors, the counterparts 
of Pr,l = 1/2(1 ± 75) for tensor currents [46J. In 
the following, we will be interested in the difference 



n 



TT 



n 



TT' 



This correlator is an order pa- 
rameter of SxSB, thus its high energies can be param- 
eterized entirely by OPE condensates. Since the chiral 
tensors F± (q) are quadratic in the momentum, the 
situation is analogous to the vector channel, namely 

km iW) = £ . (59) 



K-oo) 



n — 3 



n TT ( g 2 ) 



fl 



-q 2 + m 



(63) 



v 



The expansion of Eq. (f6"3"|) at high energies and subse- 
quent matching to the OPE of Eq. (|59|) leads to 



iti? - fl 



2A, 



fW B = , 



IB'"-B 
1 2 4 
f B m B 

2 „,6 



S8 • 



(64) 



Notice that, contrary to Ulr, there is no pion pole 
in the tensor correlator. In its place one finds the 
tensor susceptibility (3f|, which is related to A3 as 
Xt — ^^(qq)^ 1 ■ Eqs. are formally identical to 
Eqs. (fTT)) for the vector channel 8 and therefore the pre- 
diction for the low energy parameter fl is 



A = i 

32 



m 2 



1 6 



8 (Ci) 



T 

A 2 



(65) 



Unfortunately, and contrary to Liq, very little is known 
about fl, the main reason being that neither fl nor A3 
are accessible to experiment, and only lattice QCD or 
sum rule analyses can provide an estimate. Dividing 
Eq. (fT5")) by Eq. ([55)) and applying factorization, one 
gets the prediction 



-Ego 

J 71 



3 fl 
4A§ 



(66) 



which in principle could be tested in lattice simulations. 
The previous equation can be compared with the VMD 
prediction. Using the values of Ref. [38] for the decay 



7 Similar observations were made in Ref. |4;I and especially l44fl 
in the study of the same correlator with a toy model with 
QCD ultraviolet constraints. More recently, and in the same 
direction, there have also been claims that LMD does not seem 
to work for the scalar sector 14511. 



However, the nature of and A3 is very different: 
nothing to do with the breaking of chiral symmetry. 



A3 has 
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couplings, together with raj = 1.234 GeV, m p 
GeV and f B = 0.18 GeV, one finds 



Ji 3 



J VMD 



A 3 



0.770 



(67) 



£s = C£f, where C = — 8L w f^ 4 . With this relation 
one can then perform a fit for ^6 and £sj resulting in 



This disa gree ment seems to comply with the claims in 
Refs. 0, Sl| against the reliability of VMD. 



VI. CONCLUSIONS 

The use of meromorphic functions to describe cor- 
relators in QCD is nowadays a standard procedure at 
the base of methods like VMD, QCDSR or the MHA. 
From a formal standpoint, all these methods can be 
viewed as phenomenological approximations to QCD 
in the large number of colors. However, for a long time 
no clear distinction was made between Euclidean and 
Minkowski applications. This might have been the rea- 
son why there are some claims in the literature that the 
large- N c limit of QCD is unrealistic. There is no doubt 
about this if one is confined to the physical axis. How- 
ever, this point of view is not justified for Euclidean 
observables. Actually, as long as there is Euclidean 
matching, Euclidean observables are very insensitive 
to changes in the Minkowski axis. 

This phenomenon was first observed in studies of the 
MHA, where the fundamental importance of Euclidean 
matching was emphasized, but only recently the con- 
nection with Pade approximants has been made. In 
Pade theory, QCDSR and MHA are physical realiza- 
tions of partial-Pade and Pade-type approximants, and 
this stability of Euclidean quantities is a natural con- 
sequence of Pommerenke's theorem. 

Following Pade theory, the main philosophy adopted 
in this work is that meromorphic ansatze are purely in- 
terpolators, with poles and residues that need not bear 
any resemblance to physical parameters. The method 
proposed is a plain Pade approximant, where poles 
and residues are constrained in a maximal way. The 
method is Euclidean in the sense that (a) no hadronic 
input is used and (b) predictions should be restricted 
to Euclidean space. This allows to get analytical re- 
sults between low and high energy parameters. The 
main results are 



Ik 

HI 



J 7V 



Xo = 



(68) 
(69) 



£ 6 = (+7.6 ±0.4) ■ 10~ 3 GeV 6 ; 
£s = (+8.3 ± 1.0) • 10~ 3 GeV 8 



(70) 



It is worth stressing that the previous results give a 
self-consistent picture, with Am, = (4.60 + 0.28) MeV, 
C 87 = (4.0 ± 0.2) • lO- 3 GeV" 2 , to 2 , = (1.13 ± 0.06) 
GeV 2 and compatible with factorization. On the other 
hand, typical values for to 2 , in Eq. (|69[) yield values 
which favor a large xa-> m excellent agreement with the 
determinations of Refs. [27l[37T|. 

Bearing in mind that these results are just the first 
iteration of a Pade approximant to the Wlr and Hvt 
correlators, this self-consistency is tantalizing and a 
strong indication that Pade approximants to Hlr and 
ILvt might converge very fast. This convergence in any 
case is restricted to Euclidean space. We have seen that 
the residues and poles determined through matching 
can become complex, meaning that there is no possi- 
ble identification between the parameters in our ansatz 
and physical hadronic parameters. This is a limitation 
of the approach, but by no means pathological. The 
ansatz mimics the Euclidean behavior of the full spec- 
trum with a finite set of poles. Thus, not surprisingly, 
the effective poles so determined might not look like 
physical ones. This again follows from Pommerenke's 
theorem. 

Pade approximants might play a fundamental role 
in an eventual understanding of LMD. At this point, 
and from the results found in this paper, it seems rea- 
sonable to speculate that LMD is closely linked to the 
ultraviolet behavior of the correlators and so should 
be assessed on a case-by-case basis. For instance, it is 
well-known that LMD works reasonably well for Hlr- 
This might be due to the superconvergence proper- 
ties of Hlr- In contrast, for the scalar sector, where 
the convergence is decreased mainly by the absence of 
gauge invariance, there have been claims that meson 
dominance might not be a good approximation [45| . 
and actually the detailed study of Ref. [44[ on the low 
energy parameter Lg also points in that direction. One 
typically notices a seesaw effect between low and high 
energy parameters: when superconvergence is at work, 
low energies can be reliably determined, while high en- 
ergies are extremely challenging. This happens for in- 
stance with TLlr- In contrast, for n^r, where no su- 
perconvergence is at work, high energies can be reli- 
ably estimated while low energies are hard to handle. 
Therefore, the discrepancies found in the literature for 
Xo might be a consequence of misusing LMD. 



The first expression not only gives a prediction for £s to 
be positive (because Lio is negative) but, as Figure 1 
illustrates, it describes correctly all the phenomenolog- 
ical determinations consistent with a positive £g. This 
suggests a quadratic relation between the condensates, 
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APPENDIX A: DIMENSION 8 OPERATORS IN 
THE BACKGROUND FIELD METHOD 

I will outline the computation of dimension-6 and 
dimcnsion-8 operators in Ulr, Hs-p and Utt, to- 
gether with their values in the factorization approxi- 
mation. 

In the chiral limit, dimension-6 operators are either 
purely gluonic or purely fermionic, while dimension- 
8 operators can be purely gluonic or mixed. For the 
purpose of this paper one can ignore the two-quark and 
the purely gluonic operators altogether, because they 
cancel in the correlators to be considered, and one is 
therefore only left with the four-quark operators. In 
order to compute its contribution, we will work in the 
external field method formalism, an extensive review 
of which can be found in Ref . j47| • 

We start by expanding the full quark fields to first 
order in the gluonic field: 

u(x) 

d(x) 



(P + 



keeping the terms quadratic in P. The momentum 
operator P satisfies 

(x\P\y) - iD x 8{x ~ y) , (A6) 

where the derivative acts on the right, i.e., on the Dirac 
delta and any other quark fields depending on x. The 
strategy to follow therefore consists in using integra- 
tion by parts to isolate the Dirac deltas. All integra- 
tions can be performed in this way, with the derivative 
operators acting on the quark fields. 

Since we are only interested in the scalar parts of 
the correlators, we can project them with the following 
expression 

IW<7 2 ) = -\g^ R (q) , (A7) 



where I will be using the following conventions: 

G^ = ~[D^D V ] ■ 
9s 

D» = dn-ig^Ql . (A2) 

The gluon field can be split into a classical static back- 
ground field G^(x) upon which a dynamical quantum 
field <7 M (x) will propagate: G^{x) = G^{x) + g fl (x). Us- 
ing this decomposition in the QCD Lagrangian one can 
compute the quark and gluon propagators 

S(q) = J tfxe^WjrlO) = J d*x(x\-±j\0) , 

DfM = I d?'x{x\ T — T2 1 |0) . 

J (P + q) g a pb ac - 2g s f abc G° af3 

(A3) 

Inserting the previous equations in the generic correla- 
tor 



I 

for the vector sector, and 

n£ T (<? 2 ) = ^F^%)n™ pX { q ) (A8) 

for the tensor correlator. 

Finally, we need to perform an average over momen- 
tum, according to the following formula 

p2n 

(Pm ■ ■ -iW> = 2 n (n + iyS 9fllf12 ' " 5 ' t2 "-^ 2 ™ +P erm -) 

(A9) 

The previous formula takes into account that the met- 
ric tensor is symmetric, so that the allowed permuta- 
tions are those that cannot be reduced using the sym- 
metry properties of the metric. In the following we will 



(x) + ig s / d 4 x u (0) (x) S(x - x) 



d^\x)+ig a / d 4 x S(x - x) G^x)^ d^\x) , 



n ri r 2 (g)=i / d i xe lq - x {0\u{x)Y l d(x)d{Q)T 2 u{Q)\Q) , 

(A4) 

where Ti^ stand for generic Dirac matrices, we will 
(Al) expand the propagators in powers of q/P as follows 



J 



E 



i i m i i m i m i 

- - T f T + -rf-rf-r 



l = g*p?? 

q) 2 g a p8™ - 2g s f abc G b a0 q 2 



ga/3$" 



(P ■ q) 2 

4 q e 9 a f3S 



P 2 



-g a po + 



2g s f abc G a(3 
9 4 



(A5) 
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use the short-hand notation 



2 ' 



f£* try 1 x fAV 



"TT • (AW) 



For the dimension-6 contribution, the computation is 
straightforward and it reduces to one operator. The 
results for the different sectors are 



AV) 

^6 

AS) 

4 



4Tra s (uT«„d)(dT^u) , 
-16na s (uT a d){dT a u) , 
WTra s (uT a "/ 5 d) (dT a "f 5 u) 



the operators for the axial channel it suffices to make 
the replacements {u, d} — > 7s{w, d}. Similarly, the op- 
erators for the pseudoscalar sector can be readily found 
by replacing {u, d} — > i7s{w, d}. 

A standard strategy to estimate the value of the dif- 
ferent condensates resulting from the dimension-8 basis 
operators is to use the factorization hypothesis. In this 
approximation (which is also consistent with the lead- 
ing order in the 1/N C expansion), four-quark operators 
are proportional to the quark and mixed condensates. 
The master equation to be used is 



(All) ((«Tf d)(dTJu)) = -\c N ( (u O u)Ti(d <g> d)T 



where the axial and pseudoscalar sectors can be easily 
inferred by replacing {d} — ► 75 {d} and {d} — ► 175 {d} 
in and respectively. 

The dimcnsion-8 contribution is more involved. In 
the next page we list the full basis of operators for the 
different sectors, where G stands for the dual gluon 
field strength, defined as G M „ = l/2e At! yApG Ap . In ar- 
riving to the set of operators use has been made of the 

well-known expression D 2 — g s /2 cr^G^" + If) to- 
gether with the equations of motion for the quark fields 
in the chiral limit. In order to determine the associated 
Wilson coefficients, we will define £s = ^na s T)jOj, with 
the rjj coefficients given below: 



vY = - 


5 
~18 ; 




7 

_ 12 ; 


v 

% = - 


17 

18 ; 


if = 


-2i; 


v 

m = 


7 
9 : 


vi = 


1 

~2' 


V 

m = 


1 

9 : 


rfi = 


13 

_ 12 ; 


V 

% = - 


5 

~3 ; 


vl = 


11 

12 5 


v 

% = 


1 

3 : 


v! = 

rfi = 


1 

~2 ; 

-i; 





1 

~~ 4 5 


vt = 


1 

~2 5 


r h 


1 

~ 4 5 


nt = 


1 

_ 2 ; 


% 


5 

~ 2 5 


vi = 


-3; 




1 

~ 2 5 


vi = 


-3; 




i 

~ ~3 ! 


nt = 


2i 
V 




7 

~ 36 ; 


Ve = 


1 

_ 2 ; 


V7 


1 

~ 9 5 


4 = 


= 0; 


% 


1 

~ 6 5 


vt = 


1 

_ 3 ; 


% 


1 

~ ~4 ! 


4 = 


1 

3 ; 




1 

" ~i~2' 


V10 = 


1 

3 ' 

(A12) 



The results for the vector channel agree with the ones 
reported in Ref. [H, [l?], El] ) U P to changes of basis. To 
the best of my knowledge the scalar and tensor oper- 
ators had not been determined before. In order to get 



(A13) 

where I\j are generic Dirac matrices and 
Cat = 1 — N~ 2 is a color factor. In the previous 
formula the spinor fields can be generalized to include 
derivatives acting on them. The relevant spinor tensor 
products are 



(g<8>§) 
(D^D^q <g> q) 



32 3 a ^ 



(qGq) 



(A14) 



where the short-hand notation G = a llv G^ 1 ' has been 
used. Applying Eq. (|A13|) on Eqs. (|A11[) , one obtains 



AV) _ 
AS) _ 

4- 



AA) 
"6" 



4Tra s C^(uu}(dd} , 
—6ira s CN(uu}(dd) , 
-2na s C N (uu)(dd) . (A15) 



For dimension-8 operators we have to be more spe- 
cific: as shown in Ref. [l6| , application of the equations 
of motion and the factorization formulae do not com- 
mute, leading to an ambiguity of order N~ 2 . Here we 
will follow the prescription adopted in Ref. [l7| , which 
leads to 



e (V) _ AA) 
S8 — S8 



AS) 

S8 



C-p) 



4 = -Cs 



= —2ira s g s CN((uGu}(dd) 
= 6ira s g s CN((uGu)(dd} + 

= -^gsC N ((uGu)(dd) 4 



h (uu)(dGd)) , 
(uu)(dGd)) , 

(uu)(dGd)) . 

(A16) 



For all the applications in the main text, I will work in 
the isospin limit. Therefore, 



(uu) = (dd) = (qq) , 
(uGu) = (dGd) = (qGq) 



(A17) 
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a. Vector sector 



o\ 

oX 
o\ 
oX 
oX 
oX 



(uD 2 T a ^ 5 d)(dT^ 5 u) + (uT a ^ 5 D 2 d)(dT^ 5 u) + (uY a ^ 5 d)(dD 2 Y^ 5 u) + (uYfad)(dT^ 5 D 2 u) 
{uD^D v Y^ al5 d){dY^u) + (uT^idD^D^T^u) + (uT^ 5 D v d){dT v al5 u) + (uT^ 5 d)(dT^ 5 D V D '„«) 
(ulS,Y a x d)(dD^Y x u) 

(urS^y)(dD v T^u) 

(uD^Y a xl5 D»d)(dY x al5 u) + (uY a xl5 d)(dD„Y x l5 D»u) 



(uG^r;d)(dr a ul5 u) + (ur;G^d)(dr a ul5 u) + (nr« 75 rf)(JG^r» + (uY a ul5 d)(dY;G^u) 



(A18) 



b. Scalar sector 



Of = (uD 2 Y^d)(dY^u) + (uY^D 2 d)(dY^u) + (uY« u d)(dD 2 Y^u) + {uY^d){dY^D 2 u) 

Of = (uY^d)(dD»Y a D"u) + (uD»T a D v d)(dT^u) 

of = (uD x Y;J x d)(dY^u) + (uY^d)(dD x Y^D x ^ 

Of = (uD v D^ x d)(dY a vX u) + (uT^D^D v d)(dTl x u) + («T£* d){dD v Y a vX u) + {uY» x d){dT a vX D v D ^u) 

Of - (uD,D p T^d)(dTl x u) + («rf D v D„d)(dY a „ x u) + {uT^d){dD v D „rj A u) + («T^ dX^A^u) 

O e s = (oSTJ A B,,d)(^ A «) + («T^d)(dD M r^D wU ) + (uD^ x D^d)(dr: x u) + (uY^ x d)(dD^ a x D^u) 

Of = {UD a Y a d){dD a Y a u) (A19) 



c. Tensor sector 



Of = (uL» 2 r a (i)(dr a u) + (ur a D 2 d)(dr a u) + (uY a d){dD 2 Y a u) + {uY a d)(dY a D 2 u) 



O2 = 

°3 = 

Ol = 

Ol = 

Ol = 

07 = 



Or 

Os 



{uD 2 Y alb d){dY a lb u) + {uY alb D 2 d){dY a lb u) + {uY al5 d)(dD 2 Y a l5 u) + (uY al5 d){dY a l5 D 2 u) 
(uDY a Dd)(dY a u) + (uY a d)(dDY a Du) 
(uD^ a x l5 D»d)(dT« xl5 u) + (uTZ x d)(dDT a Du) 
{uD^D„d)(dT a u) + {uY a d){dD^D v u) 
(uD^D l/ d){dY a u) + {ur a d)(dD^D„u) 
{UD a T^d)(dD a T a ^u) 

(u*D a T^d)(dD p Y a pft u) + {UD p Y^d){dD a Y^u) 



O 
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{uD^D v d)(dT a ^ 5 u) + (uY al5 d)(dD„Y a ap D u u) 
g s [(uG^T^d)(dT a u) + (uY^G^d)(dY a u) + (uY a d)(dG^Y a ^u) + {uY a d){dY a ^G^ u)] 

{uG pu Y^d){diY a lb u) + {uY^G^d){diY a lb u) + {uiY a lb d){dG^u) + (mY a l5 d)(dY^G^u) 



(A20) 
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